Improving meson two-point functions by low-mode averaging 
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Some meson correlation functions have a large contribution from the low lying eigenmodes of the Dirac operator. 
The contribution of these eigenmodes can be averaged over all positions of the source. This can improve the signal 
in these channels significantly. We test the method for meson two-point functions. 



1. Introduction 

Simulations with light quarks face the problem 
of increasing noise in meson two-point functions 
with decreasing quark masses. Some of this noise 
is due to using only a single source for the meson 
instead of averaging over all possible positions of 
this source. We present a method to improve such 
meson correlation functions. More details can be 
found in 1 ; for related work see |2] and [H] . The 
idea is to average the contribution of the low- 
lying eigenmodes of the Dirac operator D over all 
positions of the source on the lattice. The contri- 
bution of the high modes is calculated using the 
standard method, inverting the Dirac operator on 
a single or a few quark sources. Ref. P] has in- 
troduced the name low mode averaging (LMA) 
for this method. LMA can have a considerable 
impact on the noise in the two-point function if it 
is dominated by the low- lying modes. We expect 
it to be particularly effective at low quark masses 
and for Dirac operators with good chiral symme- 
try like the overlap operator 0] which we use in 
our test. 

An example for the dominance of the low 
modes is the pseudo-scalar scalar difference 
shown in Fig. ^ The contribution of the low ly- 
ing modes almost saturates the signal in the re- 
gion 7 < t < 30, from which the masses are to be 
extracted. 

For the simulation we have done, the improve- 
ment comes essentially for free since the eigen- 
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Figure 1. The PP-SS two-point function. The 
contribution from the 20 lowest eigenmodes is 
shown. 



modes are used to precondition the inversion of 
the Dirac operator. The associated gain in speed 
of the inversion alone justifies their computation. 

2. Method 

We are considering zero momentum two-point 
functions of the form 

= E trriG(x',t"-t-t;a;",t")x 

x' ,x" ,t" 

T2G{x",t"-x',t" + t)) , 
with G the propagator matrix, 

{D + m)G{x\t';x,t)^5l,,,5tM , 



2 



and Fj the Dirac matrix corresponding to the me- 
son in question. 

In a standard simulation, however, one does not 
compute the fuU propagator but only a row, by 
inverting the Dirac operator on a localized source 
at, e.g., {x',t') = (0,0). On a given set of N 
configurations one thereby obtains an estimator 
(C^it)) for {C{t)). 



{Ci{t)) 



^^trriG(0,0;a;,t)x 
r2G(x,i;0,0)) 



(1) 



However, using eigenmodes of the Dirac opera- 
tor, one can compute C{t) itself. As this is in gen- 
eral not possible for all eigenmodes, we split the 
propagator in a contribution from the low modes 
and one from the high modes G = Gl + Gh with 
Gl the propagator in spectral representation 



GL{x,t-x',t')^Yl 



{x,t\j){j\x't') 
iXj + m 



(2) 



The sum is over the n lowest eigenmodes of the 
Dirac operator \j) with eigenvalues iXj. 

The two-point function is thereby separated 
into four different parts, one from the low modes 
alone, one from the high modes and two interfer- 
ence terms 

G{t) = GLL{t) + CHL{t) + CLH{t) + GHH{t) . (3) 

The first contribution Gll can be expressed by 
the low-lying eigenmodes alone and can thus be 
averaged over all positions of the source. The 
other contributions are restricted to the usual one 
or a few quark sources. 

3. Test of the method 

To test the method we use 80 quenched gauge 
configurations of size 12^ x 36 generated with 
the Wilson gauge action at /3 = 5.9. We use 
the overlap Dirac operator constructed from HYP 
smeared gauge links. The low mode averaging is 
done using the 20 lowest eigenmodes of the Dirac 
operator which are also used to precondition the 
computation of the propagator. The bare quark 
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Table 1 

Labels for Dirac structure of currents. 
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Figure 2. The effective mass plots for the 
PP~ SS correlator at iriq = 0.025. The open cir- 
cles are the effective mass without low mode av- 
eraging, the the signal given by the filled squares 
is averaged using 20 eigenmodes. 



masses are between 0.015 and 0.035 which corre- 
spond to pseudo-scalar to vector meson mass ra- 
tios m-ps/m\ ranging between about 0.4 to 0.64, 
see ID for further details. On each of the config- 
urations the inversion of the Dirac operator was 
done on two Gaussian sources with radius 3a, one 
located on time-slice i = 0, the other on i = 16. 
We average over these two positions. The F ma- 
trices depend on the meson in question and we 
shall use the abbreviations given in Tabled 

In Fig. 12 we plot the effective mass for the 
pseudo-scalar scalar difference correlator. (Sub- 
tracting the scalar two-point function removes the 
contributions of the zero modes.) Using LMA 
with 20 eigenmodes, we observe a significant im- 
provement in the noise of the signal. In Ref. ^ 
we show similar plots for other correlators. 

To give an impression of the improvements in 
the various channels, the ratio of the uncertainty 
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Figure 3. Ratio of the error bars for n — 20 com- 
pared n = at time-slice t = 5 for the different 
meson correlators. Data uses two source points 
for the high eigenmode part of the correlator. 




Figure 4. Dependence of the uncertainty in the 
extracted PP — SS mass on the number of eigen- 
modes included at anio — 0.025. 



at time-slice t = 5 with LMA over 20 eigen- 
modes to the uncertainty without LMA is shown 
in Fig. El In the 75, 757^, 707^ and 7^ channels, 
we find an improvement of roughly 30%, which 
corresponds to a factor of 2 in statistics. The 
7570 and 7i7j cannot profit from the LMA. It is 
known that low modes make a small contribution 
to these correlators 0. 

The ultimate goal of this method is to improve 
the errors of the extracted meson masses. In 
Figs. 21 and |S1 we show the dependence of the er- 
ror of the masses on the number of eigenmodes 




Figure 5. Same as Fig. ^ for the Boi mass. 



used in the LMA. We find that the improvement 
seems to saturate at 12 to 16 modes included and 
reaches about 30% for the PP-SS mass and about 
40% for the vector meson from the Boi- 

4. Conclusion 

We tested a method to improve meson two- 
point functions in lattice QCD. We found an im- 
provement in the uncertainty of some masses cor- 
responding to a factor of 2 in statistics. This im- 
provement comes at virtually no additional cost 
because the cost of computing the low eigenmodes 
is justified by the acceleration of the inversion of 
the Dirac operator alone. 
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